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Abstract 

We have investigated the dynamics of vortices in a dilute two dimensional Josephson 
junction array where a fraction of the superconducting islands is missing. We have 
used the multiple trapping model to calculate the mobility of vortices and the 
frequency dependence of the resistance and inductance of the array. 
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We consider an infinite two dimen- 
sional (2D) Josephson junction array 
(JJA), where each site is randomly 
occupied by a superconducting is- 
land with some probability p {p = 1 
for regular array) and empty with 
probability 1 — p. At a lower crit- 
ical value of p = pc {pc = 0.5 for 
triangular array, for site percola- 
tion) the superconducting proper- 
ties [1,2] of the arrays are destroyed. 
This type of arrays undergoes a 
transition from superconducting to 
normal states at a critical temper- 
ature which is called Berezinskii- 
Kosterlitz-Thouless transition tem- 
perature {Tbkt)- Experimentally [3] 
the frequency dependence of the con- 
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ductance G = 1 / Rg + 1 / iflLg of the 
sheet (array) is observed {fl is the 
frequency). We calculate the con- 
tribution to the resistance (Rg) and 
inductance (L^) produced by moving 
vortices (V) and antivortices (A). For 
this we have to first find the mobility 
yU of the V and A in the disordered 
array and for it we use the multiple 
trapping model (MTM), developed 
for electronic transport in amorphous 
semiconductors [4]. A vortex moves 
randomly in the array and whenever 
it reaches a hole (area of missing 
sites) it experiences some pinning 
potential by the hole and falls into 
it. With some probability it can also 
be released from the hole due to 
thermal excitations. We have consid- 
ered the holes of missing sites in 2D, 
for simplicity, to be circular. In this 
model the vortex mobility /i is given 
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(z/'(c<j) + iv"{uj))^Q where v{uj) = 
= i^'iuj) + iu'^uj). /io is the mo- 
bihty of V or A in a regular JJA. 
D{N) is the number density of holes 
with N missing sites taken from 
some experimental fit [5] and Nt is 
the total number of sites of the array. 
E{N) = l^2kBTBKT'^^VNTT is 
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the binding energy in a circular hole, 
j3 = T is the temperature. 
In our model we get the following 
expressions for Rg and Lg 
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where A' = ^'n^^^^jni''{uo) and 
A" = 87i^^nu"(uj), t = n is 

P-Pc t ^ ' ' Tbkt ' 

the vortex density, u = Q/ua is the 
scaled frequency and uja{^ 10^ Hz) 
is the Debye frequency of the lattice. 
-Ro and Lq are the resistance and in- 
ductance resulting from the currents 
flowing in the junction in absence of 
vortices. 

Three types of dependences of Rg on 
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Fig. 2. L-^iu)) vs a; for different T. 

p (p ~ 1) but for lower p {p pc), 
X rises to 2, respectively 4/3, for 
T = Tbkt respectively T = IOTbkt- 
For inductive part oc at low 
uj while at high u, Lj^ oc u leav- 
ing an intermediate frequency range 
where is independent of u. The 
critical frequency for crossover from 
oc to constant decreases 
with the decrease of T and increase 
of p. Our results are in good agree- 
ment with the experimental data [3]. 
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Fig. 1. Rs{uj) vs LO for different T. 

UJ are visible. At very low and high uj 
white spectrums are observed. Inbe- 
tween Rg oc uj^ with a; = 1 for higher 
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